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ABSTRACT 

Explosive hydrogen burning in type I X-ray bursts (XRBs) comprise charged 
particle reactions creating isotopes with masses up to A ~ 100. Since charged 
particle reactions in a stellar environment are very temperature sensitive, we use 
a realistic time-dependent general relativistic and self-consistent model of type I 
x-ray bursts to provide accurate values of the burst temperatures and densities. 
This allows a detailed and accurate time-dependent identification of the reaction 
flow from the surface layers through the convective region and the ignition region 
to the neutron star ocean. Using this, we determine the relative importance of 
specific nuclear reactions in the X-ray burst. 



Subject headings: X-rays: bursts — nuclear reactions, nucleosynthesis, abun- 
dances — stars: neutron 
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Introduction 



Type I X-ray Bursts (XRBs ) (see iBildstenl <ll998h: IStrohmaver fc Bildstenl (120061 ) for 
reviews) were first explained by Woosley &: TaamT ( 1976 ) who associated the XRBs with 
thermonuclear runaways on the surface of neutron stars that accrete a mixture of hydrogen 
and helium from a semi-detached low mass companion star (jJosslll977l ). 



In Woosley & Taam's thermonuclear flash model the impact of the accreting matter 
fully ionizes the matter and heats it to 1 — 2 x 10 8 K explaining the persistently observed X- 
ray emissions. The accreted matter then undergoes gradual compression as it sinks as new 
matter continuously piles on top of it. Under these atmospheric conditions the electrons 
are degenerate, but the nucleons are not. There fore the matter is subject to a thin-shell 
instability that triggers a thermonuclear runaway (IHansen fc van Hornlll975l ). For accretion 
rates below roughly one Eddington (M = 1.12 x 10 18 gs _1 ), the bottom layer of the newly 
accreted matter becomes unstable after a few hours/days and burns explosively and gives 
rise to a burst of X-rays as the atmosphere is heated to 1-1.5 GK. The sudden release of 
nuclear binding energy heats the atmosphere rapidly and increases the luminosity within a 
few seconds after which the luminosity decays exponentially as the atmosphere coo ls again 
producing the observational features of a type I X-ray burst (jJosslll978l ; lTaamlll980l ). These 
bursts are the most common thermonuclear explosions in the universe, so the indestructable 
and repeatable LMXBs are useful for indirect conclusions about the behavior of matter under 
extreme conditions. 



X- ray bursts have bee n explored theoretically by lHanawa fc Sugimotol (119 83); 



(119871 ); iKoike et al.l (119991 ) and stable burning has been explored by lSchatz et al 



Fuiimoto et al 



(11999 ) us- 



ing relatively simple hydrodynamical models to estimate the burning conditions e.g. a set 
of (p,T,X), where p is the density, T is the temperature, and X is a composition array 



describing the fractional concentration of each isotope. On t he ot 
1 dimensional models have been constructed by many groups (Joss) 



1979c lHanawa fc Sugimotol Il983l ; IWallace et al.lll982l ; lAyasli fc Joss 



ler hand more realistic 
1978) : iTaam fc Picklum 



19821 ) but they suffered 



from relatively simple nuclear reaction networks so only recently m odels have successfully 
included both aspects JWooslev et al.lbooi [Fisker et al.ll2o"oi . l2006h . 



1981); 



fll998h : 



he relevant types of reaction sequ e nces in XRBs have been d i scussed bvlWa 



Champagne &: Wiescherl (119921 ); Ivan Wormer et al.l (119941 ) ; lHerndl et al.l (119951 ) ; ISchatz et al 



lace &; Woosley 



Schatz fc Rehml (120061 ). Important are (p, 7)-, (a, 7)-, (p, a)-reaction rates as well as 



/3-decay rates between the valley of stability an d the proton driplin e. Reaction rates must 
be known up to the end-point of the rp-process ( iSchatz et al.ll2001al ). 



In the past serveral attempts have been made to identify critical reaction rates in X- 
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ray bursts. It is, however, not possible to directly test the as tronomical number o f possible 
perturbations of the thousands of participating reaction rates. I Woosley et al.l (12004 ) changed 
group s of decay rates and narro wed the rates down to several important candidates. iFisker 
(|2004h : lFiBker et al.l tool 200g) relied on "inspired guesses" and found individual important 
rates. Recently, lAmthor et al.l (120061 ) has used a one-zone model and individually varied a 
large number of reaction rates with the intent of verifying "one-zo ne" -candidates with a 
multi-zone model. Using Monte Carlo methods, iRoberts et al.l (120061 ) varied random groups 
of reaction rates and similarly identified the most significant candidates for later verification 
with multi-zone models. 

In this paper, we use such a full ID X-ray burst model with a complete nuclear reaction 
network to answer one of the fundamental questions in this field: what are the nuclear 
reactions that power X-ray bursts? While previous studies have used simplified models to 
delineate basic types of reaction sequences we can now go the next step and describe the 
actual nuclear reaction sequences that occur as a function of time and depth during a typical 
X-ray burst. Because temperature, density, and initial composition vary greatly as a function 
of depth there is no single reaction flow, but a range of very different sequences that influence 
each other. Identifying the nuclear reactions that take place in X-ray bursts is a prerequiste 
for understanding X-ray burst light curve features in terms of the underlying nuclear physics 
and for determining the nuclear physics uncertainties in X-ray burst model predictions of 
lightcurves and other observables. It is also essential to guide experimental and theoretical 
efforts to address these uncertainties in the future. 

Cross sections have typically been predicted by global models which in most cases have 
been fitted to stable nuclei and extrapolated to proton-rich nuclei. In many cases cross sec- 
tions have also been predicted by nuclear shell model calculations. However, with upgrades 
of existing expe rimental facilities and new facilities many of these reactions are now in rang e 
of experiments fcappeler etaDl 19981 : Iwiescherll200ll : Iwiescher k SchatzllioOll ISchatzlbooj ). 



With a better understanding of the nuclear physics it will also become possible to 
address potential issues beyond the ID approximation, such as the interplay of lateral flame 
propagation and nuclear energy release timescales that can affect the modeling of burst rise 
times. 



2. The ID multi-zone computational burst model 

In this paper we compute and describe one XRB model using the parameters M = 
1.4M , R = 11km and a proper global accretion rate of M — 1 ■ 10 17 g/s. This choice 
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yields a H/He-ignited XRB corresponding to case (1) of iFujimoto et al.l (ll98ll ) where the 
observation al data depends more on the rp-process than is the case for a pure He-ignited XRB 
case (2) oflFuiimoto et al.l (119811 ) ) . The Newtonian acc r etion rate (observable at infinity) of 



Schatz et al.l (12001a 



b) and model zM of IWooslev et all (120041 ) was M = 5 ■ 10 16 g s -1 , so the 



burst behavior of our model is expected to be similar but have a slightly higher hydrogen 
content and a slightly lower helium content at the point of ignition resulting in a longer rise 
time. 

This model is calculated using a general rel a tivisti c type I X-ray burst simulation code 
that is described in more detail in Fisker et al. (12006). The c ode couples the general rel- 
ativistic hydrody namics code, AGILE (iLiebendorfer et al.ll2002l ). with the nuclear reaction 
network solver of iHix fc Thielemann! (Il999i). The code includes radiative, conductive, and 
convective heat transport as described in iThornd (119771 ) and uses an arbitrarily relativistic 
and arbitrarily degenerate equation of state. We calculate the radiative opa cities due to 



homson scattering and free-free absorption using the analytic formulations of ISchatz et al. 



(119991 ) . We use the same co nductiv i ty for mulations for electron scattering on electrons, ions, 



phonons, and impurities as iBrownl (120001 ) . 



AGILE solves the general relativistic equations in a spherically symmetric geometry on 
a comoving grid. The computational domain covers about 7 pressure scale heights and is 
discretized into 129 log-ratioed grid zones with a column densitjQ ranging from y — 1.2 x 10 6 
g cm -2 to y — 3.9 x 1 9 g cm" 2 . The c omputational domain is bounded by a realistic 
core boundary interface (iBrown 2003 . 2004 ) and a relativistically corrected grey atmosphere 
(jThornelll977l ; IWeiss et al.ll2004l ) which is integrated numeric ally out to P ?ur f = 10 18 gcm~ 2 
using a 4th order Runge-Kutta method for greater accuracy (IFisker et al.ll2006l ). 



The rp-process is naturally limited once it reaches the A ~ 104 region because neutron 
deficient nuclei in this mass range become a-u nbound. This termi nates the reaction flow 
via (7, a) reactions and forms a SnSbTe cycle (ISchatz et al.l l2001al ). This determines the 
maximum network size that i s needed, unless the a lpha unbound nuclei can be circumvented 
in multiple proton exposures. ISchatz et al.l (j2001af ) demonstrated that the A ~ 104 endpoint 
can be re ached if burst peak te mperatures and hydrogen concentration at ignition are high. 
However, IWooslev et al.l ( 20041 ) showed, that these ignition conditions are only fulfilled for 
the first burst after the start of accretion on a pure 56 Fe atmosphere. Compositional inertia 
effects (ITaaml Il993l ) for subsequent bursts significantly reduce peak temperature and the 



1 The relativistic column density is mass of a column above an area: y = J R _ r p^" wnere T = 
\J\ — 2GM/Rc 2 , so P ~ gy, where R is the neutron star radius, M is the neutron star mass, p is the 
density, P is the pressure, and g — GM/TR 2 is the surface acceleration of gravity. 
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amount of hydrogen at ignition as ignition occurs at a lower pressure and depth thus limiting 
the rp-process to A < 64 with only a fraction of heavier isotopes produced. 

The nuclear reaction network used in this work employs 304 iso topes (see table []])■ A ll 
the connecting particle reactions are taken fr om the REACLIB (see ISakharuk et al.l (120061 )). 
These reaction rates have also been used in IWeinberg et al.l (120061 ) . The network includes 
all isotopes between the valley of stability and the proton dripline up to 64 Ge. Here isotopes 
with /5 + -half lives > 1 day are considered "stable" on the timescale of the burst intervals , 
so their daughters are not included. The hot proton-proton chains of IWiescher et al.l (119891 ) 
are also included. Above the 64 Ge waiting point, only isotopes between the proton drip 
line and half lives less than 1 minute are included. This is because protons only capture 
on these high-Z isotopes during the burst's peak temper ature which i s only sustai ned for 
a few seconds. Weak rates up to Z = 32 are taken from iFuller et al.l Jl98ol . Il982al lbl) and 



Langanke fc Martmez-Pinedol (120011 ). Since only a small fraction of material is processed 



above Z = 32, it is a reasonable approximation to ignore neutrino losses from heavier isotopes 
(jSchatz et al.lll999l ). These considerations significantly reduce the size of the networkm which 
decreases the simulation run-time. The inner boundary, i.e. towards the neutron star crust, 
has been slightl y impr oved compared to previ ous work which used either a massive substrate 



(I Woosley et al. 



2004)) or pa r amet er values (jRembgesI (119991 ) ) . This work uses the neutron 



star core code of lBrownl (120001 . 120031 ) which calculates the thermal luminosity emanating from 
the crust given the temperature at the atmosphere-ocean interface. The code includes pair, 
photo, and plasmon neutrino emission. This neutrino luminosity is only a few percent of 
the thermal luminosity but still several orders of magnitude larger than the hydrodynamical 
uncertainty due to the conservative formulation of mechanical equations. Different types 
of convection occur when thermal fluctuations cause instabilities to grow. Their rate of 
growth determine the eddy-velocity, so all instabilities can be treated by the mixing length 
theory (MLT) implementation. The present work only includes the Schwarzschild-Ledoux 
instability, because it is the dominant form of convection during the burst (which is the 
only period relevant to this paper), whereas secular instabilities (e.g. semi-convection) occur 
in between bursts and are negligible at high accretion rates, because the diffusion speed 
is smaller than the advection speed of the accretion. The initial model was computed by 
running the simulation for hundreds of bursts until the burst ashes had advected completely 
to the bottom of the model and the computational envelope was in thermal balance with the 
neutron star core model. At this point, the envelope was considered to be self-consistent, 
that is, independent of any possibly unphysical initial values, and a typical X-ray burst was 
picked for analysis. 
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Table 1: The table shows the list of isotopes which describes the rp-process. See the main text 
for details. An earlier version of this re ac tion n etwork has been used in the following works 



(jFiskerj 12004!: iFisker et al. 



200i l2005al fl hoOfih but it now includes the hot p roton-proton 



chains of IWiescher et al.l (119891 ). The network is described in more detail in IFisker et al. 
(l2006h . 
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Burst simulations 



Fig. [T] shows the luminosity for the analyzed burst as a function of time. When com- 
paring it to observations, it should be kept in mind that our model assumes a spherical 
symmetric ignition whereas in reality ignition most likely happens at a single point on the 
neutron star after which the flamefront spreads and eventually covers the entire neutron star. 
Therefore Fig. [I] may be thought of as the luminosity of a single point under the assump- 
tion of negligible latera l heat transport. How ever, this assumption seems to provide a good 
comparisson to reality (jGalloway et al.ll2004l ). 
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Fig. 1. — The figure shows the luminosity as seen from an observer at infinity as a function 
of time for a typical burst. The timescale has been reset so that t — corresponds to the 
peak luminosity of the burst. 
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The luminosity is the combined product of the energy released from nuclear reactions 
at different depths. The Kippenhahn diagram in Fig. [2] shows the specific nuclear energy 
release rate as a function of column density and time as well as the extent of the convective 
zone. 

The ignition region is easily indentified as a sudden and rapid increase increase in nuclear 
energy generation. This causes a temperature spike that causes a convective instability. Since 
convection transports heat very efficiently, the zones above the ignition point also ignite. 
This is clearly seen in Fig. [2J However, this heat transport quickly restores the shallower 
radiative/conductive temperature gradient and convection quickly ends. The Kippenhahn 
diagram also demonstrates how nuclear energy generation decreases as fuel depletes and how 
residual helium from the previous burst contributes to nuclear reactions below the ignition 
region. Still, most of the nuclear energy is released, not at the point of ignition but in the 
hydrogen rich layers immediately above the ignition region. 



Reaction flow 



Runaways occurring in a mixed H/He layer mainly proceed via the rp-process ([Wallace fc Woosley 



( 1198ll )). where the characteristic timescal e, r rp ~ y)Ti/9, is given by the sum of the half- 



lives of the /3-decays in the reaction flow (Ivan Wormer et al.l (119941 ) ) . However, depending 



on the flow pattern, a simultaneuosly occurring (a,p)-process, which does not depend on 



, „, — ~, 

Schatz et al. 


(1998 


); 


Fisker et al. ( 



ignited XRB, the temperature gradient only produces a minor convective instability. 

The analysis of a one-dimensional X-ray burst model is very complex, as the model is 
characterized by rapidly changing temperature conditions and nuclear reaction sequences 
in each layer. All these effects are tightly interconnected through energy generation and 
heat transport by radiation and convection. Since the different layers interact and also burn 
differently due to different compositions and temperature, the burst can not be understood 
based on the burning of one layer only, but must be analyzed for several different burning 
layers. Therefore the analysis is split into four regions: the region around the ignition point, 
the convective region, the surface, and the ocean, which are sufficiently different to merit 
separate attention. These regions are shown in Fig. [3] which shows a trace of the burst 
conditions for different depths (pressures) during a complete revolution of the limit cycle. 



Following the cooling of the previous burst, the individual layers reach their lowest 
temperature and highest density of the cycle. The subsequent accretion increases the hy- 
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Fig. 2. — The figure shows the specific nuclear energy release rate (color coded) as a function 
of logarithmic (log 10 ) relativistic column density (see footnote CD) and time as well as the 
extent of the convective zone (black sail shaped outline). The black circles correspond to the 
descriptions of the reaction flow in £JD Starting from the top of the figure, they are: ocean 
(14. 5p . ignition region (14.1 .2| - fl~l .6[ 14.1.11 is not shown), above ignition region (|4.2.2l - H~2.9p . 
14.2.11 and 14.2.101 are not shown), bottom of the convective region ( 14.3.1fl4~3.2l) . and surface 

040). 
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Fig. 3. — From left to right (solid line): y = 1.8 x 10 6 g/cm 2 (surface), y — 1.0 x 10 7 g/cm 2 
(top of the convective region), y = 2.5 x 10 7 g/cm 2 , y = 3.7 x 10 7 g/cm 2 (bottom of the 
convective region), y = 7.9 x 10 7 g/cm 2 (above ignition), y = 9.5 x 10 7 g/cm 2 (ignition point), 
and y = 1.2 x 10 8 g/cm 2 (ocean). The dashed line indicate the region which is convective 
during the rising of the burst. The circles and their associated numbers correspond to the 
figures in §H 
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drogen content of the layer, which in turn lowers the density, because the incre ased electro n 



abundanc e of hy drogen requires less mass to maintain the hydrostatic pressure (IJosd (119 77); 



oss fc Lil (119801 ) ) compared to the heavier and more neutron-rich ashes (lHanawa &; Fujimoto 
( 119841 ) ) . This is most clearly seen in Fig. [3J just above the ignition region, which decreases its 
density by about a factor two during the quiescent phase, as the electron-rich surface ashes of 
the previous burst sink into this region. For an accretion rate of M = 1 • 10 17 g/s and a recur- 
rence time of At « 11000s, the neutron star accretes a mass of AM = M ■ At » 1.1 ■ 10 21 g 
(5.5 • 1O~ 13 M ) in between bursts. This means that matter above y = 5.8 x 10 7 g/cm 2 
is freshly accreted, whereas matter below comprise the old surface ashes of the previous 
burst (s), therefore the composition in the ignition region actually consist of heavier ashes 
with a comparably lower hydrogen/ helium abundance. When the matter ignites and the 
nuclear runaway causes a rising temperature, it eventually affects the degeneracy of the elec- 
trons and decreases the density further so the trace runs up the left leg of the cycle in Fig. [3J 
until the fuel is exhausted as it burns into heavier ashes shortly after the peak temperature 
is reached. The /3 + -decays during and subsequent to the rp-process decrease the electron 
abundance and bring the trace down the right leg as the envelope cools. Therefore the sep- 
aration in density between the rising leg and the decaying leg accounts for the change in 
composition, so the largest change happens around the ignition regions, whereas the surface 
does not change its composition much. The different compositions and hydrostatic pressures 
with corresponding temperatures and densities of the burning regions change dynamically on 
a nuclear timescale defined as min(dt / d In Yi), where Yi is the abundance of the ith isotope. 
Therefore the analysis of the nuclear reaction flow proceeds in a different way compared 
to previous works, which assumed solar abundances burning a t fixed densities and tem - 



peratures and describ ed the integrated flow over many minutes ( Ivan Wormer et al.l (119941 ); 



Rembges et al.l (119971 )); instead the instant flow rate is described as the thermodynamic state 



variables change. 

The net reaction flow rate from isotope % to isotope j is defined by 

fij — ~fji = Yi^j ~ > (1) 

where Y^j is the time rate-of-change of the abundance of the ith isotope resulting from all 
reactions converting isotope i to isotope j. The flow-rates for the different times of Fig. [3J 
will be demonstrated in the flowcharts of the following sections, which describe the reaction 
flow rates at the ignition point, the region above it, the convective zone, the surface (of our 
model), and the ashes going into the ocean. 

In these figures the main reaction-flow is described by the heavy lines. Very thick 
lines indicate (p, 7)(7,p)-equilibrium. Thin lines indicate a flow rate just above 10 _6 mol/g/s 
increasing their thickness logarithmically to a maximum after which they stay constant. Also 
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shown are the mass fractions, Xa = ^A=A-^-ii f° r a gi ven mass, A, as a function of mass 
number. If Xa > 0.20 for any A, the bar is cut off and replaced with a dotted line. This 
only happens for the A = 1 and A = 4 cases where the concentration can be read in the 
figure caption. 



4.1. Ignition region 

Between bursts, the surface ashes from the previous burst sink down under the weight of 
the newly accreting matter while the hot CNO cycle transforms hydrogen into helium. The 
hot CNO cycle is beta-limited and therefore burns at a fixed rate that mainly depends on 
the concentration of 14 and 15 0. It is partially moderated by a quiescent breakout via the 
15 0(a, 7) 19 Ne reaction which depletes 15 and th us slows down the conversion of hydrogen 



into helium. This is discussed in more detail in iFisker et al.l (120061 ). As seen in figure [5j 
the reaction flow can subsequently return to the hot CNO cycle via the bi-cycle 19 Ne((3 + u) 
(T 1/2 = 17.22s) 19 F(p,a) 16 0(p, 7 ) 17 F(p, 7 ) 18 Ne(/?+,z/) (T 1/2 = 1.672s) 18 F(p,a) 15 which 
speeds up the conversion of hydrogen to helium a nd thus influences the com position for the 



runaway. This cycle is discussed in more detail in lCooper fc Narayanl (120061 ) . 



The hot CNO cycle increases the 4 He concentration until a runaway of the extremely 
temperature sensitive triple-alpha reaction ensues and causes a spike in the nuclear energy 
release rate. 



The triple-alpha reaction creates 12 C which immediately captures two protons to become 
14 0, causing the fraction of 14 (Ti/ 2 = 76.4s) to increase as seen in Fig. HJ since 14 
{Ti/2 — 76.4s) decays too slowly. Meanwhile the increasing temperature of the nascent 
nuclear runaway leads to a breakout of the hot CNO cycle into the rp-process. The details 
are described in the following subsections. 

We now describe the reaction flow in terms of temperature, density, and proton and 
alpha fractions as the time develops. The time is synchronized, so t = coincides with the 
peak surface luminosity. 



4.1.1. Fig. 0- T = 2.86 ■ 10 8 K, p = 6.81 ■ W 5 g/cm 3 , X = 0.09, Y = 0.42, t 



-103.078s 



At this t ime t he increasing temperature has caused the flow rate of 15 0(a,7) 19 Ne (see 
Fisker et al.l (120061 ) for a detailed discussion of this rate) to reach 10% of the 15 0(/5 + , i/) 15 N 
rate establishing a breakout of the hot CNO cycle (the 1% limit was breached at t = —558s) 
which extends into the light iron region. 
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Fig. 4. — The hydrogen, helium, CNO type matter, and metal (the rest) mass fractions as a 
function of time. The time scale has been synchronized to coincide with the burst luminosity 
peak at t — 0. Notice the run-up in 14 immediately prior to the runaway. Also note that 
the hydrogen in the ignition region is completely exhausted during the burst, while about 
~ 5% helium remains. 
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At this point proton captures establish a flow out of 19 Ne. The matter in this flow can no 
longer return to the hot CNO cycle and the reaction flow proceeds with 19 Ne(p, 7) 20 Na(p, 7) 
21 Mg, where it is blocked by photodisintegration, because of the 21 Mg(p, 7)(7,p) 22 Al-equilibrium. 
Therefore the flow proceeds via 21 Mg(/3+ v) {T 1/2 = 0.124s) 21 Na(jo,7) 22 Mg(/3+,z/) (T 1/2 = 
3.32s) 22 Na(p, 7 ) 23 Mg(p, 7 ) 24 Al(p, 7 ) 25 Si. 



Here the flow branches into either 25 Si(/3 + ,z/) (T 1/2 = 0.198s) 25 Al(p, 7 ) 26 Si (T 1/2 = 
1.84s) (p, 7 ) 27 P or 25 Si(p, 7 ) 26 P(p + ,u) (T 1/2 = 0.020s) 26 Si(p, 7 ) 27 P or 25 Si(p, 7 ) 26 P(p,7) 
27 S(/3 + , v) 27 P all of which have 27 P as the end point. The characteristic time depends on the 
mass fraction weighed harmonic mean of the half lives of the beta decays along respective 
pathways. As the temperature rises, the proton capture branches become initially more 
dominant but then decrease again as photodisintegration of the weakly proton bound and 
short lived proton-rich P and S isotopes steers the flow away from the dripline again. Yet 
at this particular temperature the beta decay path of 25 Si dominates the proton capture to 

26 p 

From this point, 27 P(/? + , v) (T 1/2 = 0.242s) 27 Si(p, 7) 28 P(p, 7) 29 S(/?+, v) (T 1/2 = 0.146s) 
29 P(j», 7) 30 S which has a half life of 1.07s. The Q-value of proton capture on 30 S is only 
290.6keV which makes 31 C1 subject to photodisintegration at higher tem peratures. At late r 



times, this can have a large effect on the observed luminosity of the burst (IFisker et al.ll2004j ). 



Here, the proton capture still dominates, so the flow proceeds via 30 S(p, 7) (Ti/ 2 
1.07s) 31 Cl(/3 + ,z/) (Ti /2 = 0.270s) 31 S. This (T 1/2 = 2.13s) isotope either bet a decays and 
return s to 28 Si via 31 P(p, a) 28 Si or captures a proton and proceeds via 31 S(p, 7) (jlliadis et al. 
([l999b) 32 Cl(/3+ z/) (T,n = 0.285s) 32 S(p, 7 ) 33 Cl(p, 7 ) 34 Ar(/3+,z/) (T 1/2 = 0.811s) 34 Cl(p, 7 ) 
35 Ar(p, 7 ) Jlliadis et all Jl999h ) 36 K(/5+,z/) (T 1/2 = 0.302s) 36 Ar(p, 7 ) 37 K(p, 7 ) 38 Ca. 



Since 39 Sc and 40 Sc are almost proton unbound the flow must wait for 38 Ca {Ti/ 2 = 
0.416s) and 39 Ca (Ti/2 = 0.799s) to /5 + -decay before the flow stops at the well-bound 40 Ca iso- 
tope. A CaScTi cycle exists on the well-bound 40 Ca so that 40 Ca(p, 7) 41 Sc(p, 7) 42 Ti(/5 + , v) 
(T 1/2 = 0.189s) 42 Sc(p, 7 ) 43 Ti(/5+,z/) (T 1/2 = 0.429s) 43 Sc(p,a) 40 Ca. The breakout from 
this cycle happens from 43 Sc which proceeds to capture protons going through 44 Ti and 45 V 
until it ends at 52 Fe. 



The total timescale for this sequence is (c.f. lvan Wormer et al.l (119941 )) r = ln(2) _1 Ti/ 2 
8s, which i s slower than the time it takes to cover the star with a deflagration wave by a 
factor four (IFryxell fc Woosleylll982l ). Therefore a one- dimensional approximation is still rea- 
sonable. Later when the (a,p)-process ignites and the temperature increases, the reaction 
flow will move closer to the dripline decreasing the /5-half-lives, thus making the timescales 
comparable. At that point, our model is no longer fully predictive of hydrodynamically in- 
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Fig. 5.— Ignition: T = 2.86 • 10 8 K, p = 6.81 • 10 5 g/cni 3 , x = 0.09, Y = 0.42, t = -103.078s. 
(see end of f]3]for an explanation of the diagram). 
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fluenced (extensive) observables such as the time-dependent luminosity. However, our model 
still provides a local (intensive) description of the burning conditions, and therefore a realistic 
description of the reaction flow. 



4.1.2. Fig. 1T = 3.99 ■ 10 8 K, p = 6.41 ■ 10 5 g/cm 3 , X = 0.05, Y = 0.36, t = -12.938s 

Approximately 90 seconds later the 14 0(a,p) 17 F reaction reaches 1/3 of the flow rate 
of the 15 0(a, 7) 19 Ne-reaction. This starts a second hot CNO bi-cycle: 14 0(a,p) 17 F(p, 7) 
18 Ne(/3 + ,v) 18 F(p, a) 15 which runs alongside the bi-cycle discussed above in §4.11 

At this stage 22 Mg(p, 7) 23 A1 and 22 Mg(j3 + ,v) (T 1/2 = 3.34s) 22 Na become comparable. 
Consequently the flow path through 22 Mg(p, 7 ) 23 Al(p, 7 ) 24 Si(/?+,z/) (T 1/3 = 0.190s) 24 Al 
competes with 22 Mg(/3 + ,v) (Ti/ 2 = 3.34s) 22 Na(p, 7) 23 Mg(p, 7) 24 A1 effectively creating a 
shortcut. Since the flow rates are about equal, the effective timescale becomes the flow rate 
weighted harmonic mean of the two half-lives ~ 0.10s, which is much faster than before. This 
reduces the total timescale to reach Ca to ~ 5s. A similar shortcut exists with 25 Si(p, 7) 
26 P(p, 7 ) 27 S(/3+, u) (T1/2 = 0.021s) 27 P competing with 25 Si(/5+, v) (T 1/2 = 0.188s) 25 Al(p, 7 ) 
26 Si(p, 7) 27 P however, here the proton capture Q-value is only 141keV, so the faster path is 
reduced by photodisintegration. 

At this time, the concentration of 31 CI has peaked and is now being destroyed by photo- 
disintegration. Therefore the flow must pass through the 30 S(/? + , v) (Ty 2 = 1.08s) reaction, 
which is the slowest weak reaction in the flow and adds about a second to the total timescale. 

Reaching 31 S the flow now branches again. Instead of going through the slower 31 S(/? + , v) 
(T1/2 = 2.13s) 31 P(p, 7) 32 S(p, 7) 33 C1, the flow can now go directly through either 31 S(p, 7) 
32 Cl(/3+,z/) (T 1/2 = 0.293s) 32 S(p, 7 ) 33 C1 or 31 S(p, 7 ) 32 Cl(p, 7 ) 33 Ar(/5+,z/) (T 1/2 = 0.153s) 
33 C1 which shaves another 2 seconds off the characteristic time for the rp-process. 

The flow now breaks into the p/-shell nuclei by proton-captures on 39 Ca and 40 Ca 



(IWiescher fc Gorred (119891 ) ) . The very fast /? + -decays on the highly radioactive Sc and Ti 



isotopes cause the flow to spread (Figj6j) and makes an analysis of the timescales difficult. 

The hot CNO like cycle discussed in the previous section now has proton capture break- 
outs on via 42 Ti(p, 7) 43 V and 43 Ti(p, 7) 44 V and 43 Sc(p, 7) 44 Ti. 

These reaction pass through 43 V and 44 V and after several combination of proton capture 
and beta decays, the flow goes through the 45 V bottleneck which can only happen through 
either a /3 + -decay (7i/ 2 = 0.59s) or a proton capture to 46 Cr. The next bottleneck is 48 Cr 
which can be reached from 46 V by either 46 V(/3 + , v) {T\ji = 0.429s) 46 Ti(p, 7) 47 Ti(p, 7) 48 Cr 
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Fig. 6.— Ignition: T = 3.99 • 10 8 K, p = 6.41 • 10 5 g/cni 3 , X = 0.05, Y = 0.36, t = -12.938s. 
(see end of §Hfor an explanation of the diagram). 



-18- 



or 46 V(p, 7 ) 47 Cr(/3+,z/) (T 1/2 = 0.497s) 47 V(p, 7 ) 48 Cr or 46 V(p, 7 ) 47 Cr(p, 7 ) 48 Mn(/3+, v) 
(T 1/2 = 0.030s) 48 Cr. 

The 48 Cr bottleneck has a half life of Tu^ = 2.02h which makes it "stable" on the 
timescale of the burst. The 48 Cr(jo, 7 ) 49 Mn reaction is therefore important at this stage 
because it is the only way for the flow to proceed. 

After capturing a proton the flow proceeds from 49 Mn to 50 Mn via either a beta decay 
followed by a proton capture or vice versa. The flow from 50 Mn to 51 Mn proceeds in a similar 
manner. The 51 Mn isotope has a half life of 7\/ 2 = 35.3m so 51 Mn captures a proton and 
becomes 52 Fe. 

There is a slight flow out of 52 Fe that moves to 56 Fe via a series of proton captures 
followed by beta decays. At this point 56 Fe captures several protons to 59 Cu which decays 
and proton captures and decays to 60 Ni. There are also proton captures on Ni-Zn ashes from 
the previous burst at this time. 



4.1.3. Fig. 0" T = 5.37 ■ 10 8 K, p = 5.81 ■ 10 5 g/cm 3 , X = 0.03, Y = 0.31, t 



-10.631s 



At this point the 14 0(a,p) 17 F rea ction is 5 times stronger t h an the 15 0(a , 7 ) 19 N e 
reaction. This st arts the (a,p)-proces ([Wallace fe Woosleyl (1198 IT ); ISchatz et al.l (119981 ); 
Fisker et al.1 (120041 )) which here runs as u O(a,p) 17 F(p, 7 ) 18 Ne(a,p) 21 Na(p, 7 ) 22 Mg(a,p) 



2". 



Al(p,7) 26 Si. 

Another (a,p)-reaction exists on 21 Mg(a,p) 24 A1. This reaction soon overpowers the 
22 Al /3 + -decay from the 21 Mg(p, 7 )( 7 ,p) 22 A1 equilibrium which becomes largely irrelevant 
for the burst from this point. 

While 31 C1 above 30 S is in (p, 7 )( 7 ,p)-equilibrium with 32 Ar, the main flow goes through 
the beta decay of 30 S which has a significant impact on the energy generation as it blocks 
the rest of the flow of therp-process. 

In the Ca-Fe region, the flow moves closer to the dripline with 43 V(p, 7 ) 44 Cr(/3 + ,i/) 
(T 1/2 = 0.030s) 44 V and 47 Mn(p, 7 ) 48 Fe(/?+, u) (T 1/2 = 0.030s) 48 Mn becoming more active. 

Above Fe, the flow extends to 60 Zn as more proton rich nuclei start to capture protons. 
With 60 Zn's Ti /2 = 4.10m halflife, 59 Cu's T 1/2 = 83.1s halflife, and the stable 58 Ni, the 
60 Zn(p, 7 ) 61 Ga reaction is the only way to move the flow forward. 



In addition, proton captures on heavier isotopes from the previous burst moves towards 
the dripline. This compositional inertia increases the average mass and charge of the final 
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Fig. 7.— Ignition: T = 5.37 • 10 8 K, p = 5.81 • 10 5 g/cni 3 , X = 0.03, Y = 0.31, t = -10.631s. 
(see end of §Hfor an explanation of the diagram). 
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ashes. 

4.1.4. Fig. 0- T = 7.30 -10 8 K, p = 5.07 ■ 10 5 g/ cm 3 , X = 3.5 x 10" 4 , Y = 0.26, t = -9.980s 

The protons are now almost exhausted. This can also be seen in Fig. HI Once this 
happens, the proton-rich isotopes near the driplines decay towards the valey of stability, 
where they undergo (a,p)-reactions. The figure shows the situation with less than 1% 
hydrogen so proton captures are still occurring although they are weakening. Also, hydrogen 
from (a,p)-reactions can still serve as a catalyst in the (a,p)-process. 

The (a,p)-process extends further and now includes the 25 Si(a,p) 28 P and 26 Si(a,p) 29 P 
reactions. 

SiPS, PSC1, ArKCa, and CaScTi cycles are observed with 28 Si, 32 P, 36 Ar and 40 Ca as 
the nexus, but they are not consequential to the flow. 

In the Zn region, the flow proceeds via 59 Cu(/3 + , v) (Ti/ 2 = 84.5s) as the proton captures 
on 60 Zn remain weak. This flow reaches 63 Ge and 64 Ge but will not move further in this 
region as the hydrogen concentration is rapidly depleting. 

4.1.5. Fig. 0- T = 9.03 ■ 10 8 K, p = 4.46 ■ 10 5 g/cm 3 , X = 3.3 x 10" 6 , Y = 0.12, t = -8.075s 

At this point the region receives more heat from adjacent regions than it produces. This 
allows endothermic (a,p)-reactions on better bound nuclei. 

Presently, sufficient material has been moved to the 60 Zn isotope ensuring its decay 
and extending the rp-process from light isotopes into heavier isotopes. However, the only 
protons available for capture on these heavy nuclei come from (a,p)-reactions on A = 20-36 
so the general lack of protons ensures that the rp-process at this point does not proceed 
beyond A = 64. Additionally the shortage of protons means that the flow moves away from 
the proton dripline with the remaining protons generally capturing on the currently most 
abundant nuclei (now determined by half-life) with the largest cross sections and the lowest 
Coulomb barriers, that is, A = 20-36. 
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Fig. 8.— Ignition: T = 7.30 • 10 8 K, p = 5.07 • 10 5 g/cm 3 , X = 3.5 x 10" 4 , K = 0.26, 
t = —9.980s. (see end of $1] for an explanation of the diagram). 
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Fig. 9.— Ignition: T = 9.03 • 10 8 K, p = 4.46 • 10 5 g/cni 3 , X3.3 x 10" 6 , Y = 0.12, t = -8.075s. 
(see end of §4] for an explanation of the diagram). 
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4.1.6. Fig.UM T 



9.89 ■ 10 8 K, p = 
X28 



4.20 • IOV cm 3 , X 
= 0.22, t = -3.018s 



2.4 x 10~ 9 , Y = 0.05, 



The flow through the alpha-chain nuclei is clearly seen in Fig. [10] which shows the 
reaction flow at the time where maximum temperature is achieved. Notice that 12 C(p, 7) 
13 N(a,p) 16 is much stronger than the direct 12 C (a, 7) 16 Q-reactio n as long as the (ce,p)- 
reactions are still possible on heavier isotopes ([Weinberg et al.l 120061 ) . 



The reaction flow continues with (a, 7)-reactions up to 36 Ar. Eventually the downward 
heat flux from the upper regions becomes too weak to sustain the (a,p)-reactions and the 
reactions die out leaving only radioactive isotopes, which slowly decay to stabler ones. 



4.2. Above the ignition region 

It is relevant to know the reaction flow and its energy release at the depth, which reach 
the highest temperature during the burst, because it heats up adjacent and colder regions. 
This is because heat is transported as electrons and photons are diffused along a negative 
temperature gradient. The highest temperature of a burst ignited by mixed hydrogen and 
helium is reached just above the point of point of ignition. 



4.2.1. Fig. [22t T = 2.80 ■ 10 s K, p = 4.29 ■ 10 5 g/cm 3 , X = 0.41, Y = 0.43, t = -103.078s 

The 15 0(a, 7) 19 Ne- reaction is less important at this depth, because less "hot CNO mate- 
rial" has been created due to lower operating temperatures and densities of the triple-alpha 
process. So while the reaction burns off the existing 15 0, the runaway at this depth oc- 
curs, when the heat from the ignition point below increases the 3a reaction, so the 12 C(p, 7) 
13 N(p, 7) 1A 0(a,p) 17 F reaction path dominates the eventual runaway. This is because 14 
{T\/2 = 70.6s) does not have the time to decay during the runaway. However, at this point 
the 15 0(a, 7) 19 Ne-reaction does establish a very weak flow to the iron-region along a flow 
path identical to the initial path in the deeper region. 

Although this region contains former surface ashes there are no further captures beyond 
these isotopes at the present time. 
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Fig. 10.— Ignition: T = 9.89 • 10 8 K, p = 4.20 • 10 5 g/cni 3 , x = 2.4 x 10~ 9 , Y = 0.05, 
X28 = 0.22, t = —3.018s. (see end of §2]for an explanation of the diagram). 
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Fig. 11 — Above ignition: T = 2.80 ■ 10 8 K, p = 4.29 • 10 5 g/cm 3 , X = 0.41, Y = 0.43, 
t = —103.078s. (see end of J|]for an explanation of the diagram). 
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4-2.2. Fig. EE' T = 3.97 ■ 10 8 K, p = 3.68 ■ 10 5 g/cm 3 , X = 0.400, Y = 0.406, t = -12.387s 

At this point 14 0(a,p) 17 F reaction reaches 1/3 of the flow rate of the 15 0(a, 7) 19 Ne- 
reaction. This happens at a lower temperature for this depth, because the 14 0/ 15 abundance- 
ratio is relatively higher. 

The breakout establishes the second hot CNO bi-cycle (discussed in §4.1.21) ; in contrast 
to the ignition point, the first bi-cycle (discussed in §4.1 ft is not established, because it is 
already sufficiently hot and there is sufficient hydrogen to capture on 19 Ne destroying it 
immediately. 

The heat flux building up rapidly from the ignition point below means that short cuts 
e.g. 22 Mg(p, 7) 23 Al(p, 7) 24 Si(/3+, u) (T 1/2 = 0.190s) 24 Al competing with 22 Mg(/5+, v) (T 1/2 = 
3.37s) 22 Na(p, 7) 23 Mg(p, 7) 24 Al quickly becomes dominated by the leg closest to the proton 
dripline. 

From this point the reactions are identical to the flow described in §4.1.21 Since there 
is more hydrogen in this region the reaction on isotopes heaver than Mn are faster, yet since 
the temperature at this point in time (same as §4.1.2[) is lower, the the capture rates on 
lighter isotopes than Mn are slower. 

4.2.3. Fig. \£kT = 4.44 • lO 8 ^, p = 3.46 • 10 5 g/cm 3 , X = 0.398, Y = 0.402, t = -11.091s 

14 0(a,p) 17 F is now as strong as 15 0(a, 7) 19 Ne. The flow through the 31 Cl(/? + , v) (T 1/2 = 
0.270s) waiting point is currently approximately equal to the flow through the 30 S (Ty 2 = 
1.08s) waiting point, but the latter will quickly become dominant as rising temperatures 
prevent the formation of 31 CI due photodisintegration. 

In the Zn region, the flow stops at the long lived 59 Cu (Tu 2 = 84.7s) and 60 Zn (Ty 2 = 
4.29m), but the temperature is not yet sufficiently high for proton captures to established a 
reaction flow to heaver nuclei, nor has sufficient time passed to allow a substantial amount of 
material to decay through these two nuclei. It is interesting to note that processing to heaver 
material either depends on the temperature becoming sufficiently high for the 60 Zn(p, 7) (7,2?) 
61 Ga equilibrium to allow (p, 7)-reactions on 61 Ga or the temperature remaining sufficiently 
low for the flow to decay through the faster 59 Cu(/? + , v) 59 Ni-reaction. 
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Fig. 12.— Above ignition: T = 3.97 • 10 8 K, p = 3.68 • 10 5 g/cm 3 , X = 0.400, Y = 0.406, 
t = —12.387s. (see end of §4] for an explanation of the diagram). 
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Fig. 13.— Above ignition: T = 4.44 • 10 8 K, p = 3.46 • 10 5 g/cm 3 , X = 0.398, Y = 0.402, 
t = —11.091s. (see end of §4] for an explanation of the diagram). 
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4.2.4- Fig. W\T = 5.75 ■ 10 8 K, p = 2.19 ■ 10 5 g/cm 3 , X = 0.393, Y = 0.387, t = -10.418s 



At this point the 18 Ne(a,p) 21 Na-reaction activates making it possible to move into the 
rp-process via 12 C(p, 7) 13 N(p, 7) 14 0(a,p) 17 F(p, 7) 18 Ne(a,p) 21 Na and so forth instead 
of "waiting for the Tu 2 = 1.67s beta-decay of 18 Ne. The 18 Ne(a,p) 21 Na-reaction is thus 
especially important as most of the energy release in the atmosphere originate in the re- 
process at lower depths/higher hydrogen concentrations as seen in Fig. [2j 

However, presently 90% of the the flow through the lighter isotopes stops at the 30 S 
{T1/2 — 1.09s) waiting point with only a small flux following from its decay. This causes a 
temporary dip in the energy production, though the higher temperature ensures a flow close 
to the dripline from Ca to Ni. This flow is, however, slowed down at the N = 28 isotones 
due to the long half-lives of 55 Co (T 1/2 = 10.3h) and 56 Ni (T 1/2 = 24.9h) which effectively 
prevents any beta decays of these isotopes. It is also interesting to note the (p, a)-reactions 
on the heavier Cu isotopes going back to Ni while releasing hydrogen. 

At this point there is a weak flow out of 59 Cu which allows additional proton captures 
viz. 59 Ni(p, 7) 60 Cu(p, 7) 61 Zn(jo, 7) 62 Ga which can either decay or capture an additional 
proton to 63 Ge which then decays. If 62 Ga decays it captures two additional protons and 
goes to 64 Ge. It is also possible to reach 65 Ge if 63 Ga decays. 

Heavier isotopes are generated through proton captures on ashes from the previous 
burst. None of these are beta decaying though. 



4.2.5. Fig. T = 6.97 ■ 10 8 #, p = 2.51 • W^/cm 3 , X = 0.381, Y = 0.372, t = -9.994s 

A couple of seconds after its breakout, 14 0(a,p) 17 F, becomes so fast that any 14 
is immediately destroyed. Consequently 15 is only created via the hot CNO bi-cycle. 
However, the bi-cycle will become void, because it is now sufficiently hot for alpha-particles 
to penetrate the Coulomb barrier of 18 Ne, thus skipping its Ty% = 1.67s /3 + -decay. 

Additional (a,p)-captures now happen on 21 Mg and 25 Si. The latter circumvents the 
(T1/2 = 0.176s) half-life of 25 Si, thus shortening the characteristic reaction flow timescale 
in the A = 20-30 region. The timescale is dominated by the 30 S(/? + , v) (Tin = 1.09s) 30 P 
reaction since 98% of the flow passes through this reaction. 

Another reaction is 24 Si(a,p) 27 P but that is not as significant since 24 Si gets destroyed 
by photodisintegration. It is interesting to note a weak but present 18 Ne(a,7) 22 Mg which 
competes with the (a,p)-process. 
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Fig. 14.— Above ignition: T = 5.75 • 10 8 K, p = 2.19 • 10 5 g/cm 3 , X = 0.393, Y = 0.387, 
t = —10.418s. (see end of §4] for an explanation of the diagram). 
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Fig. 15.— Above ignition: T = 6.97 • 10 8 K, p = 2.51 • 10 5 g/cni 3 , X = 0.381, Y = 0.372, 
t = —9.994s. (see end of §2]for an explanation of the diagram). 
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The flow up to 58 Ni remains the same. The increased temperature and flow sets up 
NiCuZn, ZnGaGe and GeAsSe cycles on 58 Ni, 60 Zn, and 66 Ge. 

We note th at there is still no fl ow out of 64 Ge(p, 7) 65 As as 65 As is weakly proton bound 
and 2p capture ( Schatz et al. 19981 ) is not effective. Therefore the reaction flow proceeds via 
the slow beta decay 64 Ge(/3+, u) (T 1/2 = 84.9s) 64 Ga or the lighter 63 Ga(/3+, v) (T 1/2 = 26.6s) 
63 Zn. 

Reaching the N = 33 isotones, the flow reaches 67 Se (Ti/2 = 0.060s) and 68 Se (Ty 2 = 
35.5s). Further progress either depends on another 2p-reaction (jSchatz et al.lll998l ) or 68 Se 
or 67 As decaying. 



As similar challenge is posed by Kr (Tiy 2 = 17.2s), Sr (Ti/2 = 8.9s) and Zr (Ti/ 2 = 
3.9s). Presently the flow has not moved farther although protons have started captures on 
heavier isotopes of the ashes of the previous burst. 



4.2.6. Fig. EE" T = 8.34 • 10 S K, p = 2.15 • 10 5 g/ cm 3 , X = 0.358, Y = 0.346, t = -9.097s 

At this time the temperature is sufficiently high for photodisintegration of 27 S to prevent 
the shortcut, which was previously established between 25 Si and 27 P. 

However, at the same temperature the 21 Mg(a,p) 24 Al and the 22 Mg(a,p) 25 A1 reactions 
become significant. In addition, 24 Si(a,p) 27 P, 25 Si(a,p) 28 P, and 26 Si(a,p) 29 P become 
significant. Circumventing the 30 S(/3 + ,z/) (Ti/ 2 = 1.09s) 30 P becomes possible through the 
29 S(a,p) 32 Cl-reaction which at this point is not as strong as the beta decay. 

38 Ca starts photodisintegration, but since 39 Sc is proton unbound, the flow must await 
a (T1/2 = 0.04s)-decay, since this reaction is a bottleneck. 

Heavier isotopes with iV > 32 continue to 84 Mo (Ty 2 = 3.6s). Some hydrogen is 
capturing on iV = 43 and iV = 44 isotones of the ashes from the previous burst. 



4.2.7. Fig. \ThT = 8.96 ■ 10 8 ^ p = 2.07 ■ 10 5 g/cm 3 , X = 0.327, Y = 0.326, t = -8.075s 



At this time, the 28 S(a,p) 31 C1, 29 S(a,p) 32 C1, and 30 S(a,p) 33 C1 are all active. The latter 
now competes directl y with the 30 S beta decay. This competition is especially important at 
lower accretion rates (IFisker et al.l 120041 ) . 



Heavier isotopes continue to 88 Ru (T1/2 
heavier isotopes. 



Lis). There are no significant captures on 
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Fig. 16.— Above ignition: T = 8.34 • 10 8 K, p = 2.15 • 10 5 g/cni 3 , X = 0.358, Y = 0.346, 
t = —9.097s. (see end of §4] for an explanation of the diagram). 
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Fig. 17.— Above ignition: T = 8.96 • 10 8 K, p = 2.07 • 10 5 g/cm 3 , X = 0.327, Y = 0.326, 
t = —8.075s. (see end of §2]for an explanation of the diagram). 
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4-2.8. Fig. HE T = 9.93 ■ lO 8 ^, p = 2.43 ■ 10 5 g/ cm 3 , X = 0.143, Y = 0.234, t = -3.013s 

This region has now reached its maximum temperature. The reaction flow-path is very 
similar to the flow in Fig. [171 

One notable difference is the 34 Ar(a,p) 37 K reaction which is the last among the (a,p)- 
reactions for the temperatures encountered in type I XRBs. 

Additionally, heavier isotopes continue to 92 Pd and 93 Pd. This effectiv e ly cons titutes 
the end of the rp-process which is short of the prediction of ISchatz et al.l (l2001al ). The 
reason is the thermal and compositional inertia as well as the much lower peak temperature 
ac hieved by our model. If these are ignored, the flow does reach the SnSbTe cycle as shown 
bv lWooslev et al.l hoo4 ). 

This flow structure is maintained until hydrogen runs out. 



4.2.9. Fig. [M T = 9.62 • 10 8 #, p = 3.54 ■ W 5 g/cm 3 , X = 2.5 x 10" 5 , Y = 0.175, 

Xq = 0.346, t = 1.476s 

The last protons capture on the currently most abundant nuclei, namely the isotopes in 
the Ca-Ge region as the flow falls back towards the valley of stability as it is decaying along 
constant mass numbers. 

With helium still burning, the (a,p)-process is still active along with previously men- 
tioned (a,p)reaction up to A = 36. However, the 12 C(p, 7) 13 N(a,p) 16 path prevents much 
formation of 17 F. This means that the (a,p)-process starts on 21 Na which is reached from 
16 0(a,7) 20 Ne(p,7) 21 Na. 



4.2.10. Fig. \M T = 6.97 • lO 8 ^, p = 5.10 • 10V cm 3 , X = 1.0 x 10 -11 , Y = 0.089, 

X 60 = 0.347, t = 28.788s 

Half a minute after the burst the rp-process no longer operating and the temperature 
has decreased 30% from the maximum temperature. 

An (a,7)-chain connects 16 to 32 S which eliminates most of the 12 C via 12 C(p, 7) 
13 N(a,p) 16 0. Meanwhile, heavier isotopes follow constant A decay-chains back to the valley 
of stability. 
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Fig. 18.— Above ignition: T = 9.93 • 10 8 K, p = 2 A3 ■ 10 5 g/cm 3 , X = 0.143, Y = 0.234, 
t = —3.013s. (see end of §4] for an explanation of the diagram). 
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Fig. 19.— Above ignition: T = 9.62 • 10 8 K, p = 3.54 • 10 5 g/ cm3 > X = 2.5x 10~ 5 , Y = 0.175, 
X m = 0.346, t = 1.476s. (see end of §4] for an explanation of the diagram). 
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Fig. 20.— Above ignition: T = 6.97 • 10 8 K, p = 5.10 • 10 5 g/cm 3 , X = 1.0 x 10 -11 , Y = 0.089, 
Xq = 0.347, t = 28.788s. (see end of £jl]for an explanation of the diagram). 
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4.3. Convective region 

The size of the convective region is shown in Fig. [3], which shows a trace of the burst 
conditions for different depths (pressures) during a complete revolution of the limit cycle. 
Note that the convective zone only exist during the phase where the temperature rises (the 
cycle revolves clockwise). The figure shows that the convective zone does not reach the top 
of our model for this burst, but stays in a narrow region between y = 5.7 x 10 6 g/cm 2 and 
y = 5.2 x 10 7 g/cm 2 . Additionally, this burst does not reach super-Eddington luminosities, 
so no ashes will be ej ected by a radiativel y driven wind; something which is possible in 



helium-ignited bursts ([Weinberg et al.ll2006l ). 



The quantitative analysis of the turbulent convective burning is complicated by the 
mixing of matter between convective zones, which occurs as soon as and as long as a slightly 
superadiabatic temperature gradient is established. However, the convective timescale, 
T con = A/v e dd. ~ 10~ 6 -10~ 5 s <C r rp , is generally faster than the typical timescale of the 
rp-process, so the explosive burning will have alm ost the same com position throughout the 



entire convective zone (see the convective model of iRembgesl (119991 ) which assumes identical 
composition throughout the convective zone for comparison) although burning happens at 
different temperatures and densities at the top and bottom of the convective zone respec- 
tively. Furthermore, turbulent convective burning does not happen above temperatures of 
7 x 10 8 K, so the (a,p)-process, which has a much shorter timescale, does not become ac- 
tive. Therefore this region could be computed by models with a simplified description of 
the compositional evolution but a more complex (2D) hydrodynamical implementation. The 
following analysis of this region concentrates on the bottom of the convective region, because 
it is hotter and denser, and therefore the reactions proceed faster here. 



4.3.1. Fig. \M T = 3.92 ■ 10 8 K, p = 1.161 • 10 4 g/cm 3 , X = 0.610, Y = 0.363, t = -10.229s 

The convective period during this burst lasts about 2.2 s during which fresh unburned 
matter from the colder top of the convective zone mixes into the warmer bottom and back 
again. This means that temperature dependent particle-captures are effectively weaker, 
whereas the weak decays remain unaltered. 

Due to the short duration, the region attains a maximum temperature of ~ 0.7 • 10 9 KK 
while it is convecting (~ 0.9 • 10 9 K at the peak), so the high temperature short cuts like 
(a,p)-reactions or the upper leg of rp-process bifurcations available to the deeper layers and 
described in the previous sections never come into play; instead the initial reaction flow is 
mostly represented by the flow chart of Fig. [211 
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Fig. 21.— Convection region: T = 3.92 • 10 8 K, p = 1.161 • 10 4 g/cm 3 , x = 0.610, Y = 0.363, 
t = —10.229s. (see end of f]3]for an explanation of the diagram). 
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At this point, the second hot CNO bi-cycle, 14 0(a,p) 17 F(p,7) 18 Ne(/?+,z/) (T 1/2 = 
1.67s) 18 F(p,a) 15 already dominates the first hot CNO bi-cycle, 19 Ne(/3+,z/) (T 1/2 = 17.2s) 
19 F(p,a) 16 0(p,7) 17 F(p, 7 ) 18 Ne(/5+,z/) (T 1/2 = 1.67s) 18 F(p,a) 15 0, as the temperature is 
high enough for the 19 Ne(p, 7) 20 Na-reaction to dominate the slow 19 Ne(/9 + , v) {T\/% = 17.2s) 
19 F. From this point 20 Na captures another proton so 20 Na(jo, 7) 21 Mg(/3 + , v) (T1/2 = 0.124s) 
21 Na(p,7) 22 Mg. 

Here the flow bifurcates to either 22 Mg(/3+, v) (T 1/2 = 3.46s) 22 Na(p,7) 23 Mg(p,7) 24 Al 
or 22 Mg(p,7) 23 Al(p,7) 24 Si(/3 + ,z/) (T 1/2 = 0.191s) 24 Al. A similar bifurcation exists at 
25 Si(p, 7 ) 26 P(p,7) 27 S(/3+,z/) (T 1/2 = 0.021s) 27 P competing with 25 Si(/3+,z/) (T 1/2 = 0.189s) 
25 Al(p, 7 ) 26 Si(p, 7 ) 27 P. 

The 30 S waiting point still acts as a bottleneck with a small leak via 30 S(p, 7) 31 Cl(/3 + , v) 
(T1/2 = 0.272s) 30 S. However, this is quickly reduced by photodisintegration as before. 

This flow passes through the 34 Ar(/5 + ,u) (TW2 = 0.824s) bottleneck and on through 
37 Ca(/?+,i/) (T 1/2 = 0.155s) 37 K(p, 7 ) 38 Ca(/5+,z/) (T 1/2 = 0.423s) 38 K(p, 7 ) 39 Ca(/?+,i/) 
(T 1/2 = 0.808s) 39 K(j»,7) 40 Ca into the p/-shell isotopes. 

It is interesting to notice that the rp-process is already active in the 40 Ca- 52 Fe region. 
This is because the double-magic 40 Ca is the natural end-point of the reactions during the 
minute long run up to the burst ignition where a small reaction flow is already present. 



4.3.2. Fig. \M T = 6.90 • 10 8 ^ p = 9.35 ■ 10 4 g/ cm 3 , X = 0.646, Y = 0.328, t = -8.600s 

Fig. [22] shows the maximum temperature during the convective phase occuring concur- 
rently with the maximum temperature in the ignition region, which drives the superadia- 
batic temperature gradient responsible for the convective turnover. Higher temperatures are 
reached at this depth, but after this time, the region is no longer convective. 

At the end of the convective phase, the (a,p)-process at the bottom of the convective 
region extends to 26 Si. Following that, the flow to heaver isotopes is impacted by the 30 S 
{T1/2 — l-10s) waiting point. This is significant because the entire convective phase only 
lasts two half lives of 30 S. Similarly, there are bottlenecks at 59 Cu (TW2 = 91.9s) and 60 Zn 
{T1/2 — 5.28m) which require a 60 Zn(p, 7) 61 Ga breakout that does not happen at these 
temperatures during the short convection phase. 
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Fig. 22.— Convection region: T = 6.90 ■ 10 8 K, p = 9.35 • 10 4 g/cm 3 , X = 0.646, Y = 0.328, 
t = —8.600s. (see end of $4] for an explanation of the diagram). 



-43- 



4.4. Surface region 



In H/He-ignited XRBs the convective region does not extend to the top of our model. 
This means that if the convective model does not severely underestimate the convective 
strength then heavier ashes are not brought to the surface, Since the matter at the top 
of our model is extremely opaque with mean free photon paths of ~ 10~ 4 cm, the photons 
are in local thermal equilibrium (LTE) and exhibit a black body spectrum with no lines. 
Comparison between the results of this section with observations th e refore require this model 
to be coupled with a radiative transport code (see IWeinberg et al.l (120061 )). 



44.1. Fig. [H T = 5.31 • 10 8 K, p = 8.75 • 10 3 g/cm 3 , X = 0.697, Y = 0.281, t 



-0.209s 



The extent of the reaction flow at the maximum temperature is shown in Fig.[2n]and ends 
at 56 Ni. This region is limited by T < 5.3 ■ 10 8 K and the initial reactions are characterized 
by proton capt ures on the accreted heavy elements, which may have been destroyed by the 
surface impact (IBildsten et al.l 1 1 9921 ). 



4.5. Ocean (ashes) 

The inner parts of the neutron star acts as a buffer absorbing heat from the burst. 
However, for this accret ion rate it is radiated outwards again after the burst, therefore it 



does not heat the crust (IFujimoto et al.l (J1984J)). 



The early reaction flow which is caused by conductive heating in a hydrogen depleted 
environment is similar to the reaction flow in Fig. [TUJ Later it is characterized by residual 
helium, which has been advected down from the previous burst, capturing on alpha-chain 
nuclei extending to 36 Ar as shown in Fig. [2H Note that here we also have 12 C(p, 7) 13 N(a,p) 
16 being much stronger than the direct 12 C(a,7) 16 0-reaction. Here the protons are sup- 
plied by many weak (a,p)-reactions on stable isotopes resulting from long lived /? + -decays 
in the sulfur region in matter that has advected downwards from above. 



5. Conclusion 

Important in all the regions are the hot-CNO cycle and its respective breakout reactions, 
the (a, p)-process, as well as (p, 7) (7, p)-equilibria and waiting points of the rp-process. These 
are now discussed. 
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Fig. 23.— Surface: T = 5.31 • 10 8 K, p = 8.75 • 10 3 g/cm 3 , X = 0.697, Y = 0.281, t = -0.209s. 
(see end of S]l]for an explanation of the diagram). 




Fig. 24.— Ocean: T = 9.55 • 10 8 K, p = 6.06 • 10 5 g/cm 3 , X = 5.3 • 10' 11 , Y = 0.036, 
Xqs = 0.22, t = —0.209s. (see end of £@]for an explanation of the diagram). 
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5.1. Hot CNO cycles 



There are essentially three hot CNO cycles, specifically, the first hot CNO cycle: 12 C(p, 7) 
13 N(p, 7 ) 14 0(/3+,z/) 14 N(p, 7 ) 15 0(/3+,z/) 15 N(p,a) 12 C, the second hot CNO cycle: 14 0(a,p) 



17i 



F(p,j) 18 Ne(/5+,z/) 18 F(p,a) 15 and third hot CNO cycle 15 0(a, 7 ) 19 Ne(/3 



v 



19 



F(p,a) 



16 



0(p,7) 17 F(p, 7 ) 18 Ne(/3 H 



1/ 



is 



F(p,a) 15 0. 



In order to activate the second hot CNO cycle, via 14 0(a,p) 17 F and its breakout via 
17 F(p, 7) 18 Ne(a,p) 21 Na, the third hot CNO cycle must activate hundreds of second prior to 
the runaway and achieve a breakout via 19 Ne(>, 7) 20 Na. If the 15 0(a, 7) 19 Ne-reaction is too 
weak, the third cycle never activates which means that the second cycle does not activate 
either and the thermonuclear runaway does not happen (IFisker et al.ll2006l ). These rates are 
therefore quite significant in connecting the hot CNO cycle and the rp-process. Additi onally, 
the reaction flow, in particular the third hot CNO cycle of ICooper &: Narayanl (120061 ). prior 
to the runaway is important for the ignition composition as it influences the concentration 
of hydrogen and helium which is important to the thermonuclear instability. 



5.2. The (a,p)-process 



The (a,p)-process is important because it is a temperature dependent process unlike the 
rp-process that contains temperature-independent /3 + -decays. The (a,p)-process therefore 
influences the characteristic timescale of the reaction flow up to A = 36 afte r whic h the 
Coulomb barrier becomes prohibitive. Furthermore, as shown in IFisker et al.l (120041 ). the 
(a,p)-reactions in the (a,p)-process lie on waiting points with 30 S being the most signifiance. 
Other potential waiting points are 34 Ar and 26 Si. 

The most important (a,p)-reactions for the XRB are therefore 26 Si(a,p) 29 P, 30 S(a,p) 
33 C1, and 34 Ar(a,p) 37 K. The 22 Mg(a;,p) 25 Al-reaction is most likely not as important, since 
the flow moves through the 22 Mg waiting point via 22 Mg(p, 7) 23 Al before the (a,p)-reaction 
becomes active. 

Other (a,p)-reactions are less dominant since they operate at higher temperatures and 
on more proton-rich nuclei which are more susceptible to photodisintegration viz. 21 Mg(a,p) 
24 Al, 24 Si(a,p) 27 P, 25 Si{a,p) 28 P, 28 S{a,p) 31 C1, and 29 S{a,p) 32 C1. The final (a,p)-reaction 
is 13 N(a,p) 16 which in the even t of hy drogen depletion is stronger than the 12 C(a,7) 16 



reaction (also see IWeinberg et al.l (120061 )). 
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5.3. The rp-process 



The rp-process evolution depends on the concentration of hydrogen and the peak tem- 
perature. The peak temperature is easily estimated as P = a ra dT A which assumes that 
the pressure is fully supported by radiation and that the dynamical pressure is negligible. 
This is a good assumption as the gravitational binding energy is a factor ~ 20-50 higher 
than nuclear energy release of the burst. This dependence means that if the pressure of 
the ignition point is inaccurately determined, the peak temperature may be off by 10% or 
more which will significantly change the conclusions about the flow. Thermal and compo- 
sitional inerti a must be take n into account when considering the reaction flow. This was 
first shown by lWoosley et al.l (120041 ) who started wi th a pure 56 Fe a tmosp here which allowed 
accreted matter to reach deeper layers. As a result IWoosley et al. ( 20041 ) obtained the same 
results as the one- zone model of ISchatz et al.l (j2001al ) who based their ignition pressure and 
composition on analy tical estimates. On the other hand, selfconsistently obtained bursts by 



Wo os 



1999 



ev et al.l (|2004r) match t he resu lts obtained by other selfconsistent models (IRembges 



Fisker et alil2003l . I2004L l2005aE ) as well as this paper. 



5.3.1. rp-process waiting points 

Waiting points are isotopes from which further net reaction flows are (possibly temporar- 
ily) restricted due to either insufficiently high temperatures, insufficient capture particles, 
or the immediate photodisintegration due to a (p, 7)(7,p)-equilibrium. Waiting points are 
easily identified by their temporary abundance spikes. If a substantial, say 20% or more, 
part of the flow is backed up at a given isotope for a time comparable to the time scale 
of the XRB, it can significantly influence the shape of the observed luminosity curve (see 



Fisker et all ((200J)). 

During the early build up to the XRB and during the early phases of the 21 Mg(p, 7) (7, p) 22 Al- 
equilibrium, which depends on the Q-value of the proton capture reaction, means that 21 Mk 
must /? + -decay. The half life is short compared to the build-up phase which is on the or- 
der of hundreds of seconds. It is however comparable with the runaway time of the XRB. 
Therefore, the runaway depends on the 21 Mg(a,p) 24 Al reaction. Similar waiting points can 
be found along the (a,p)-process reaction path. They are 22 Mg, 26 Si, 30 S, Ar, and 38 Ca. 
The dominant waiting point in this sequence depends on the extent of the (a,p)-process 
which depends on the peak temperature. If the peak temperature is extremely high e.g. 
Tp e ak > 1-3 x 10 9 K these waiting points are bypassed by the (a,p)-process. For lower peak 
temperatures, these waiting point along with their associated Q-values and proton capture 
rates become important. However, our model has never reached peak temperatures above 
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~ 1.3 x 1 9 K for accretion rates gr eater than M — 5 • 10 16 gs -1 while accreting a solar com- 
position (lAnders & Grevessd Il989l ) on a self-consistently attained atmosphere. The 38 Ca 
waiting point might be circumvented by 38 Ca(p, 7) 39 Sc(p, 7) 40 Ti or 38 Ca(2p, 7) 40 Ti. 

Hot CNO-like cycles exist on well-bound isotopes such as 40 Ca. This isotope is par- 
ticularly interesting since the flow passes through it during the quiescent phase. The low- 
temperature 43 Sc(jo, 7) 44 Ti-reaction is therefore an important bottleneck as it determines 
the developing composition during the quiescent phase and thus the ignition conditions. The 
next bottleneck in the quiescent flow is 48 Cr(p, 7) 49 Mn. During the burst (above T ~ 5-10 8 K) 
the flow through the Ca-Ni region goes through many /3 + -decays and (p, 7)-reactions leaving 
no single determining reaction. Due to conservative scheme of numerical discretization of 
the model, it is possible to track minor variations in the luminosity due to individual rates. 

The Ni-Se region provides several waiting points. The first waiting points are 59 Cu and 
60 Zn. These half lives are on the order of the burst decay timescale and must be surpassed 
by proton captures in order for heavier isotopes to be produced. The latter is in (p, 7X7,2?)- 
equilibrium and thus depends on the Q-value of 60 Zn(p, 7) 61 Ga. There is a possible flow via 
61 Ga(jo, 7) 62 Ge. A similar situation ex ists on 64 Ge. Here 65 As is proton-unbound, so further 
flow depends on either a 2p-capture (jSchatz et al.lll998l ) or a slow /3 + -decay. This is the 
reason why most of the flow stops at the A = 64. The atmosphere cools before a substantial 
amount of matter can decay and be processed to heavier isotopes. Similar situations exist 
on 68 Se, 72 Kr, 76 Sr, and 80 Zr where the corresponding 69 Br, 73 Rb, 77 Y, and 81 Nb are also 



pr oton-unbound. These waiting points have also been identified by lSchatz et al.l (119981 ) and 
by IWoosley et al.l (120041 ) , who showed the significance of these decays by varying groups of 
electron capture and /5 + -decay rates up and down by 1 order of magnitude thus testing the 
impact of the efficiency of the reaction flow progression through the waiting points on the 
burst lightcurve. 



5.4. Superbursts and convection 



We showed that the peak burst temperature is less than ~ 1.3GK and thus not as high 
as previously assumed for nuclear reaction studies. This means that Te is not generated i n 
quantity which corro borates previous multi-zone simulations of Fis ker et al.l ( 2003. 2005b); 



Wooslev et al. 



Woosley et al 



(2004). The average mass of the ashes is ~ 64 (IFisker et al 



2003 



2005b 



20041 ) . At the same time carbon is slowly destroyed by helium captures below 



the ig nition zone and at the top of the ocean. This corroborates the findings of IWoosley et al. 



( 2004T) and it does not favor the parameter space requirements of current superburst theories 
( ICumming Bildstenl l200ll ). However, our model did not consider sedimentation effects 
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which may change this conclusion (jPeng et al.ll2006l ). 



We find that the convective region does not hit the top of our model for mixed hydro- 
gen/helium (sub-Eddington) bursters. Therefore we predict that any spectral lines observed 
during such bursts are not from material that was burned at any significant depth. How- 
ever, at lower accretion ra t es, the convective region does hit the top of our model for helium 



bursters (see iFisker et al.l (j2005al )). 



5.5. Summary 

The main result of our calculations is the identification of the nuclear reaction sequences 
that power type I X-ray bursts. In particular, we describe the complete nuclear reaction flow 
as a function of time and depth, including branchings and waiting points, as it evolves with 
realistic, rapidly changing temperatures and densities. Clearly, the reaction sequences are 
more complex than previously assumed based on the analysis of much simpler models. Our 
work is a necessary first step towards identifying the critical nuclear reaction rates in X-ray 
bursts that have the largest impact on observables such as light curves, or, indirectly, the 
composition of the ashes. One can the n also begin to disentangle the effects of nuclear burning 



on the luminosity (IFisker et al : ll2004l) from geometric effects such as the propagation of the 



burning front around the NS (jSpitkovsky et al.ll2002l ) to better explain the many different 



and somewhat inconsistent shapes of the burst luminosity profiles. 

JLF and HS were supported by NSF-PFC grant PHY02-16783 through the Joint In- 
stitute of Nuclear Astrophysics^]. FKT and JLF acknowledge support from the Swiss NSF 
grant 20-068031.02. 
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